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Analytic potentials and vibrational energies for Li 2 states dissociating to 
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Analytic potentials are built for all four states of Lb dissociating to Li(2S) + Li(3P): 

36(3^11^), 3-B(3^nu), 8 ( 7 ( 3 ^ 113 ), and 3d(3®ng). These potentials include the effect of spin-orbit 
coupling for large internuclear distances, and include state of the art long-range constants. This is 
the first successful demonstration of fully analytic diatomic potentials that capture features that are 
usually considered too difficult to capture without a point-wise potential, such as multiple minima, 
and shelves. Vibrational energies for each potential are presented for the isotopologues ®’®Li 2 , ®’^Li 2 , 
^’^Lb, and the elusive ‘halo nucleonic molecule’ ^^’^^Li 2 . These energies are claimed to be accurate 
enough for new high-precision experimental setups such as the one presented in [Sebastian et al. 
Phys. Rev. A, 90 , 033417 (2014)] to measure and assign energy levels of these electronic states, 
all of which have not yet been explored in the long-range region. Measuring energies in the long- 
range region of these electronic states may be significant for studying the ab initio vs experiment 
discrepancy discussed in [Tang et al. Phys. Rev. A, 84, 052502 (2014)] for the C 3 long-range 
constant of Lithium, which has significance for improving the SI definition of the second. 

PACS numbers: 02.60.Ed , 31.50.Bc , 82.80.-d , 31.15.ac, 33.20.-t, , 82.90.-|-j, 97, , 98.38.-j , 95.30.Ky 


Very little is known about the Li 2 electronic states dis¬ 
sociating to the 2S -I- 3P asymptote. Out of the first 5 
asymptotes (from lowest to highest: 2S -\- 2S, 2S -h 2P, 
2S'-|-3S', 2P+2P and 25-1-3^), the 25'-|-3P is the only one 
for which an empirical dissociation energy has not been 
determined for any of the electronic states dissociating to 
it. Furthermore, the only measurements that have been 
made for electronic states dissociating to 2S -h 3P, were 
for the 3c(3^E+) state [Il-@ the 67f(6^E+) state [1, 
and the 3(i(3^ng) state @-0. No measurements have 
been done on the other states dissociating to 2S -\- 3P. 

Very recently, a promising experiment has been setup 
with the ability to use photoassociation in a magneto¬ 
optical trap to make ultra-cold ®Li 2 molecules dissociat¬ 
ing to the 2S -\- 3P asymptote 0 , much like slightly ear¬ 
lier experiments which have already been successful for 
creating ultra-cold ®Li 2 molecules dissociating to 2S'-|-2P 
with very similar techniques pin. Measurements of 
the binding energies for levels very close to the 2S + 3P 
asymptote would allow for an empirical determination of 
the long-range constant which is the leading in¬ 

teraction constant in the potential energy between Li(2S') 
and Li(3P). 

At the lower asymptote of 2S -j- 2P, there is a discrep¬ 
ancy between experiment and theory for the long-range 
constant , despite Li only having 3e“ and the 

experimental value being the most precisely determined 
oscillator strength ever determined for a molecule, by an 
order of magnitude [l^. This has various consequences, 
reaching as far as limiting progress towards improving 
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the precision of the SI definition of the second [TJ . More 
precise atomic clocks are needed for various applications. 
The current definition of the second is based on a clock 
transition frequency in Cs with a relative uncertainty of 
~ 5 X 10“^®, and a commonly quoted target for improved 
precision is 10“^® [T^. The largest source of uncertainty 
limiting atomic clock precision is the blackbody radiation 
shift, which depends on the static dipole polarizability of 
the system being used for the atomic clock 0- Lithium 
is expected to play a major role in polarizability metrol¬ 
ogy, since polarizability ratios can be measured much 
more precisely than individual polarizabilities pisl] and Li 
is the preferred choice for the standard in the denomina¬ 
tor of such a ratio 0 . The discrepancy in C 3 limits the 
accuracy of a potential Li-based standard for polarizabil¬ 
ities [l3| , and hence indirectly impacts progress towards 
improving the SI definition of the second. 

Regarding the empirical value for C'g'®^^'^, for most 
electronic states, the mixing of various states towards 
the 2S + 2P asymptote significantl y comp licates the ex¬ 
pressions from which C 3 is fitted [la - l^ . The compli¬ 
cated expressions for this mixing are the same at the 
25' -I- 3P asymptote as they are for the 25 -I- 2P asymp¬ 
tote [2lj , but the fine structure splitting parameter which 
governs the significance of this mixing, is about 3.5 times 
smaller at the 25 + 3P asymptote than at the 25 + 2P 
asymptote. For 25 -|- 2P the fine structure splitting 
parameter for ®Li is AE 22 ~ ^2 ~ Pi = 
0.335 324 6 cm-i [HI, [11 [H, il|, while for 25 -f 3P it 
is only AE 32 P <_ 32 p = 0.096 cm“^ [H. Therefore, 

■j/2 -*-/2 

fj 2 S-i- 3 P jjgjgjgi; ]gg better benchmarL for an ab initio vs 
experiment comparison than , as the effect of this 

complication is smaller. 
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Measuring and assigning molecular energy levels us¬ 
ing photoassociation requires reasonably accurate pre¬ 
dictions which come from eigenvalues of a Schrodinger 
equation, and hence require a reasonably accurate po¬ 
tential energy surface. Due to the shortage of measure¬ 
ments on the Li 2 states dissociating to 2S'-|-3P, 

the most accurate potentials come from purely ab initio 
calculations. For the 3c(3^E+) and 324(3^ S+ ) states, ah 
initio calculations were reported in 1985 [^, 1995 
2006 [23 3'iid 2014 [2pJ^; and for the 3(i(3^ng) state in 
1995 U and 2014 [28l. |29||. But for the rest of the states 
dissociating to 2S + 3P; namely 36(3^nu), 6X(6^E+), 
3P(3^n„), 317 ( 3 ^ 113 ), and 6a(6^S+); the only ab initio 
calculations reported were in [^[^[^. All of these ab 
initio papers also reported potentials for states dissoci¬ 
ating to lower asymptotes, where plenty of experimental 
data is available to gauge the quality of the calculations. 

In my very recent paper on comparing experiment to 
ab initio for the 6(l^n„) state [l^ . it was found that the 
ab initio potential of [28| predicted all vibrational bind¬ 
ing energies with a disagreement of < 12 cm“^ with the 
corresponding energies of the empirical potential. Fur¬ 
thermore there was always <0.8 cm“^ disagreement be¬ 
tween the empirical and ab initio vibrational energy spac- 
ings. Finally, when comparing the dissociation energies 
J)e from to the corresponding experimental values for 
all states which have empirical S)e values available, the 
ab initio values from were never in disagreement by 
more than 68 cm“^. Therefore, the ab initio potentials 
from [23 for the states dissociating to 2S + 3P are ex¬ 
pected to be a good starting point for predicting energy 


levels with the precision required for photoassociation ex¬ 
periments as in IllL Il2l | and as may be preformed with 
the new setup in which is capable of detecting states 
dissociating to 2S -I- 3P. 

However, the ab initio calculations of [13 still have 
some major drawbacks (including, but not limited to): 

1. the ab initio points are not on a dense enough mesh 
to use as the mesh for solving the effective radial 
Schrodinger equation for predicting the vibrational 
energies {especially for large distances where the 
energies become more important for fitting an em¬ 
pirical C 3 value, and for experiments such as those 
potentially resulting from a setup such as in 0); 

2 . the ab initio points neglect the effect of spin- 
orbit coupling, which is particularly important for 
nS + n'P states of Li 2 , where the effect of interstate 
coupling has been shown to be absolutely obliga¬ 
tory for descri bing the hig h vibrational energy mea¬ 
surements [iiI[12[[17L[i3- 

3. the ab initio points do not go beyond the Born- 
Oppenheimer approximation (so do not distinguish 
between different isotopologues such as ®’®Li 2 , 
®7Li2, ^’’^Li 2 , and the elusive ‘halo nucleonic’ iso¬ 
topologues containing ^^Li), and they are non- 
relativistic. 

Drawback (1) is usually treated by fitting an interpolant 
through the ab initio points, but the resulting energies 
predicted after solving the Schrodinger equation, will be 
very sensitive to the type of interpolant used, especially 
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at the level of precision of photoassociation experiments 
(the precision in the Li 2 measurements of [ll|, [l^ was 
±0.00002 cm“^ or ±600 kHz). Also, if for sake of ease, 
a spline interpolant is used, it would be defined piece- 
wise and would have discontinuous first derivatives. The 
spline also knows nothing about the physics of nature, 
and will therefore not know what to do in regions where 
fewer ab initio points are available (in this example, for 
large internuclear distances). 


Rather than interpolating with a spline, we can fit to 
a fully analytic model potential that has the correct the¬ 
oretical behavior incorporated in the long-range region 
where fewer ab initio points are available, and this ad¬ 
dresses drawback ( 2 ) as well, since the effect of spin-orbit 
coupling at long-range can easily be incorporated into 
the model. Part of drawback (3) can also be addressed 
by fitting to a model potential, because the model can 
also build in some types of relativistic effects such as 
QED retardation, as was attempted in [Tl|, [HIE [3. 
In 2011 the Morse/long-range (MLR) potential was fitted 
to spectroscopic data for the c(l^S+) state of Li 2 where 
there was a gap of > 5000 cm“^ between data near the 
bottom of the potential’s well, and data at the very top 
3- In 2013 it was found by experiment that vibrational 
energies predicted from this MLR potential in the very 
middle of this gap were correct to about 1 cm“^ [Hj- 
Therefore, fitting the ah initio data to the MLR model 
can provide reliable energy predictions in regions where 
ab initio points are lacking or are poor in quality. 


Therefore, in this paper MLR models that incorpo¬ 
rate the long-range theoretical effect of spin-orbit cou¬ 
pling are fitted to the ab initio points from for the 
states of Li 2 dissociating to 2S ± 3P. Draw¬ 
back (3) is not addressed in this paper. However, Born- 
Oppenheimer breakdown (BOB) corrections could have 
been added to the ah initio points u sing the molecular 
electron wavefunction as described in |16l | . Alternatively, 
the entire ah initio calculation can be redone using a 
non-Born-Oppenheimer approach as has been done for 
up to 6 e“ [30|, but the a posteriori approach of doing 
a Born-Oppenheimer calculation and then adding BOB 
corrections afterwards has been shown to work better ac¬ 
cording to th e ag reement between experiment and theory 
for BeH [3l|, Also, DKH (Douglas-Kroll-Hess) rela¬ 
tivistic corrections can be added to the ab initio potential 
as in 3 1 QED effects can also be added as was done 
for H 2 in [35 and HeH+ in [^ . If any of these answers 
to drawback (3) were to be addressed by adding correc¬ 
tions to the ab initio points of [3 j fhe procedure applied 
in the present paper for fitting an MLR function to ab 
initio points, could be repeated for even more accurate 
analytic potentials. 


I. EXTENDING THE AB INITIO 
CALCULATIONS 


The ab initio calculations in did not go beyond 
22 A. However, beyond a certain length, analytic expres¬ 
sions for the potential can be derived from the theory 
of atom-atom interactions with disregard for the effect 
of overlap between each atom’s electronic wavefunction. 
These analytic expressions are based on long-range con¬ 
stants that come from atomic ab initio calculations rather 
than molecular ones, so for Li 2 the calculations only in¬ 
volve 3e“ rather than 6 e“. This means, for example, that 
a coupled cluster calculation taking account of the full 
configuration interaction (FCI) of a basis set only needs 
to go up to triple excitations (CCSDT, whose scaling 
with respect to the number of basis functions A is ~ A® 
and has been implemented since 1987 3); whereas a 
molecular calculation on Li 2 would require all the way 
up to hexuple excitations (CCSDTQPH, which scales as 
~ A^^, and has been implemented in only very few stud¬ 
ies since 2000 (STUi^ with basis sets that have not yet 
gone beyond the cc-pVDZ-DK basis set HI). Further¬ 
more, 3e“ is the limit at which the integrals have been ex¬ 
pressed analytically for explicitly correlated Slater wave- 
functions, so treating 6 e“ would either require numeri¬ 
cally calculating the integrals (which would be too slow 
even for small basis sets), or explicitly correlated Gaus¬ 
sian wavefunctions (which do not necessarily have the 
correct short- and long-range behavior). Therefore, be¬ 
yond a certain distance the analytic expressions ignoring 
wavefunction overlap but using long-range constants for 
Li based on 3e“ ab initio calculations, are expected to 
be more accurate than the 6 e“ ab initio calculations of 
3 that include wavefunction overlap. The distance at 
which this trade-off begins to lean in favor of the an¬ 
alytic expressions is heuristically given by the Le Roy 
radius 

Another advantage of using the analytic expressions, is 
that the ab initio calculations of 3 l do not include the 
effect of spin-orbit coupling, but for alkali atoms disso¬ 
ciating to nS ± n'P asymptotes, the effect of spin-orbit 
coupling at long-range has been determined analytically 
[3 3- Although all papers discussing these analytic ex¬ 
pressions to date only mention nS ± nP asymptotes, the 
expressions are also the same for nS ± n'P asymptotes 
when n ^ n' 3 l- 


A. LeRoy radii 


The 
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m-dependent Le Roy radius is given by Ji et al. 


PhR-m = 2’\/3 (^{nlm\z‘^\nlmy^^ ± {n'l'm'\z‘^\n'l'm'y^^'^ , 

_ ( 1 ) 

where for hydrogen-like atoms we have 3 l' 
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Figure 1. Point-wise original, and analytic MLR potentials for 3fo(3^nu) representing the ab initio calculations of [2^. The 
inset shows the long-range behavior in Le Roy space: it demonstrates that the original ab initio points unphysically dip below 
the theoretical curve, while the MLR behaves correctly. 



{nlrn\z^\nlrnY^'^ = 



2 3m2-Z(; + l) \ 

3 (2/-h3)(2/- 1 )) 


and for / = 0 we have (because m is also 0) 


V2 



( 2 ) 


Table I. m-dependent Le Roy radii for Li 2 electronic states 
that approach 2s -I- ml in Hartree atomic units and spectro¬ 
scopic units, and the constituent quantities that are used to 
calculate these radii. 

<r^) {nlm\z^\nlm)^t'^ rLR_m(2s + m/) rLR_m(2s-I-mi) 

a.u. a.u. a.u. A 


{nlm\z'^\nlmy^'^ = (3) 

v3 

This means that if both atoms of a diatomic molecule are 
in S states, the fact that I = I' = 0 reduces Eq. [T]to the 
original Le Roy radius of [H, : 


rLR^2((ri)V^ + (ri)V^). 


(4) 


However for a hydrogenic atom with I 0 we have 


2s 17.47 [4^ 17.47 [4^ 16.7189 8.8473 

2p 27.06 0.3810 27.0 14.0 

3p 168.69 [4^ 0.5910 40.0 20.1 


~ = 1 a.u., and using (r^) values 

calculated in [^, we are able to calculate Eq. [T] for 
Li 2 molecules dissociating to various asymptotes, and we 
present these in Table ID 


{nl\r^\niy^^ 




Ki + i)-y3 \^y^ ^ 
(5) 


where Z is the effective nuclear charge, is 

the Bohr radius scaled by the ratio of the mass of the 
nucleus mjv to the reduced mass p, of the atom. And 
for alkali atoms, the principal quantum number n is re¬ 
placed by n — a{l), where a(l) is the quantum defect 
and can be found in standard references such as Ref 
[29] of [i^. Using a(p) = 1.59 for Li, assuming that 


B. Long-range theory 

It is well-known that for large internuclear distances, 
the MLR model becomes, [13: 


V (r) ~ S)e — u{r) -!-•••. (6) 

Therefore, we can define u{r) to be the analytic expres¬ 
sion describing the theoretical interaction between the 
constituent atoms of the molecule. Each state 
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Figure 2. Point-wise original, and analytic MLR potentials for 3-B(3^ni„) representing the ab initio calculations of [2^. The 
inset shows the long-range behavior in Le Roy space. 
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considered in this paper has i^u/g daughter states result¬ 
ing from the spin-orbit coupling that lifts the degeneracy. 
For the 2u and 2g states, which are daughters of 36(3^ng) 
and 3d(3^ng) respectively, no other ft = 2 state with 


the same u/g symmetry approaches the 2S' -|- 3P asymp¬ 
tote, so the potential energy curves at long-range are not 
strongly influenced by other electronic states. Therefore, 
these states have the simplest form for u{r): 


(r) 


/ 

AE- 

V 



771=3,6,8, 

9,10,11,... / 


^3 _ ^ _ 

^3 .^6 


^8 _ ^ _ ^10 

^8 ^9 J .10 


*-"11 

,,11 


(7) 

( 8 ) 


where the zero of energy is the Li(2S'i/2) + asymptote, and AE = AE^ 2 p^^^^^ 2 p^^^ is included since the 2^ 

and 2g states both dissociate to Li(2S'i/2) + Li(3P3/2). The 36(3^nu) state additionally has a daughter state of 1„ 
symmetry, along with 3i?(3*^nii); and the 3(i(3^ng) state additionally has a daughter state of Ig symmetry, along with 
3 ( 7 ( 3 ^ 113 ). Since all l„/g states approaching 25'-|-3P have two other l^/g states of the same “/g symmetry approaching 
25' -I- 3P, the u{r) for these states is defined as the highest, middle, or lowest energy eigenvalue of the following 3x3 
matrix (including the prefactor of —1) depending on whether the state in question is the lowest, middle, or highest in 
energy respectively: 


u^’'/s(r) 


(9) 
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Figure 3. Point-wise original, and analytic MLR potentials for 3C(3^ng) representing the ab initio calculations of [2^. The 
top inset shows the long-range behavior in Le Roy space, and the bottom inset shows that the MLR successfully captures the 
tiny second minimum which has a depth of « 8.5 cm“^. 
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( 10 ) 

The notation means (^A„ or ^Ag). The zero of energy is once again the Li( 2 S'i/ 2 ) + Li( 3 Pi/ 2 ) asymptote. 

Finally, the 3b state additionally has 0+ and 0“ daughter states, and the 3d state additionally has 0+ and 0“ daughter 
states. Since all 0^ states approaching 2S + 3P have one other 0^ state of the same “/g symmetry and the same 
+/- symmetry approaching 2S + 3P, the u{r) for these states is defined as the higher, or lower energy eigenvalue of 
the following 2 x 2 matrix (including the prefactor of — 1 ) depending on whether the state in question is lower, or 
higher in energy respectively: 



/i V- cr<K2A"-'-> 

3 


m=3,6,8 

9,10,11,... 

# E 

m=3,6,8 

9,10,11,... 


C, 




3 r3 

m=3,6,8 
9,10,11,... 

-ae+i J 2 -—^ - 


m=3,6,8 

9,10,11,... 


( 11 ) 


/ 


The zero of energy is once again the Li( 2 S'i/ 2 ) + Li( 3 d^i/ 2 ) 

/ \2 

Since the leading term not shown in Eq. 0 is , 
the contribution of the C 3 terms to the long-range form 
of the potential, will interfere with the desired Cq and 


asymptote. 


Cg terms, and all Cg and Cn terms will therefore have 
spurious contributions from the cross-terms formed by 
the products of the C 3 terms with the Cq and Cg terms 
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Figure 4. Point-wise original, and analytic MLR potentials for 3d{3^Ilg) representing the ab initio calculations of [2^. The top 
inset shows the long-range behavior in Le Roy space, and the bottom inset shows that the MLR successfully captures the tiny 
second minimum which has a depth of « 13.5 cm“^. 



respectively. We fix this in the same way as was done for 
Ce and Cg in by applying a transfor¬ 

mation to all Cg, Cg, and this time also Cn terms: 

Cg ^ Cg + ^ (12) 

tig ^ Cg + (13) 

where the transformation in Eq. [T^] has to be made first 
due to Eq. [T3Is dependence on Cg. 

Additionally, the long-range formulas in terms of Cm 
constants in Eqs. ISIlOllll were derived under the assump¬ 
tion that two free atoms are interacting with each other, 
and there is no overlap of the electrons’ wavefunctions as 
would be in a bound molecule. To take into account the 
effect of electron overlap, we use the damping function 
form from : 


Cm ^ CmD^^^ (r) (15) 



where for interacting atoms A and B, p = pab = , 

in which px = y /C) is defined in terms of the ioniza¬ 
tion potentials of atom X, denoted (/^), and hydrogen 


(/^). We use s = — 1, which as shown in [i^, means that 
the MLR potential has the physically desired behavior 
V{r) ~ i/r’* in the limit as r —>■ 0. For s = —1, the sys¬ 
tem independent parameters take the values 6^1) = 3.30, 
and = 0.423 [H]. 


C. Long-range constants 

For electronic states of Li 2 approaching the 2S + 2P 
asymptote, the Cg^g^g constants for all electronic state 
symmetries have been calculated with finite-mass correc¬ 
tions for ®Li and ^Li [i^, and even an attempt at rela¬ 
tivistic corrections has been made for the Cg constants 
mill. Furthermore, for 2S + 2P, third-order pertur¬ 
bation theory has been used to calculate non-relativistic 
infinite-mass values for Cg and Cn [T^ , meaning that it 
was possible to also include the non-relativistic infinite- 
mass value of Cio calculated in . 

The situation is much less convenient for 2S + 3P. 
No third-order perturbation theory calculation has been 
done for Cg or Cn, and without Cg it does not make 
sense to include the Cig value, which was calculated in 
the same study as for the 2S + 2P asymptote [i^. Also, 
no finite-mass or relativistic corrections have been calcu¬ 
lated for the Cg^g^s values associated with 2S -I- 3P. Nev¬ 
ertheless, we have available the non-relativistic infinite- 
mass values for Cg^g^g that were calculated in [i^, and 
these were reported with an order of magnitude higher 
precision than in the very highly cited 1995 paper of 













































Table II. The best currently available long-range constants for 
Li 2 electronic states that dissociate to 2S + 3P (in Hartree 
atomic units). All values come from and were calcu¬ 
lated without relativistic corrections, and under the assump¬ 
tion that both Li nuclei have infinite mass. 




3/1 V , 



C3 

0.0033314 

-0.0033314 

-0.0016657 

0.0016657 

Ce 

3.8236x10"^ 

3.8236x10“ 

2.0282x10“ 2.0282x10“ 

Cs 

2.4870x10'^ 

2.3183x10^ 

7.8976x10® 3.7222x10® 


Marinescu and Dalgarno , and only one order of mag¬ 
nitude lower precision than the 2S + 2P values which are 
known (see Table 2 in [l^ for a list of the best known 
Cm constants for each symmetry approaching 2S + 2P). 
All Cm constants that are used in this study for 2S + 3P 
are given in Table HIl 


II. MLR POTENTIALS 


It has been suggested that fully analytic potentials , 
and specihcally the MLR may not have the flexibility 
required to capture some features such as multiple min¬ 
ima and shelves (see examples of these features appearing 
in 2S — 3P potentials of Li 2 in Figs[I]ll|). While no at¬ 
tempt (as far as I am aware) has thus far been made to 
use a fully analytic potential to capture such features, an 
increasing number of applications of the MLR potential 
after the publications of [5^[5J has made it a strong case 
for a “universal” potential form. MLR-type potentials 
have successfully described spectroscopic data for many 
electronic states of many diatomic molecules [Tl|, [Ij, 

[l^. [^. [ 4 ^ [Hhl - lT^ . It has also become customary to fit ab 
initio data for diatomic [tsI - It^ and polyatomic [H, Isol - 
systems to MLR models . Therefore, we use the MLR 
model in this study, and the results here support the idea 
of the MLR model being a strong candidate for a “uni¬ 
versal” model for potential energy curves and surfaces. 

All MLR potentials were made by fitting to the ab ini¬ 
tio points of Ref. [ 2 ^ with the program DPotFIT [s^. 
Since this is a non-linear least-squares fitting, ‘starting 
parameters’ are needed in order to allow DPotFit to 
achieve reasonable fits. Starting parameters were ob¬ 
tained from the program BetaFIT [^. When fitting to 
ab initio points, DPotFIT aims to minimize the dimen¬ 
sionless root mean square deviation: 


dd = 


\ 


-^data 

1 ^ 

^ab initioij'} 

V 

.^data ^ 


) 


where I4 ilr(*) and Vai initio{i) are the values of the re¬ 
spective potentials at the i**' internuclear distance value 
(the order of course does not matter) and A^data is the 


total number of ah initio points to which the MLR po¬ 
tential is being fitted. Uat iniUoi}) is the uncertainty in 
the i*'^ ab initio point, so that the MLR potential is likely 
to lie more closely to ab initio points at distances where 
the ab initio calculation is expected to be more reliable, 
and the requirement for the MLR potential to match the 
ab initio is less harsh in areas where the ab initio calcu¬ 
lation is expected to be less accurate. 

It is extremely hard to determine accurate estimates 
on the uncertainties for ab initio points. The ab initio 
points we are using from [ 2 ^ were all calculated with 
the same basis set (which the authors denoted by ANO- 
RCC-F), therefore there is no indication of the size of 
the basis set error. Furthermore, all of their calculations 
were done with the same number of excitations included 
in their coupled cluster method: FS-CCSD(2,0) only in¬ 
cludes I- and 2-electron excitations, so it would be ex¬ 
tremely unlikely to estimate the deviation from the full 
6-electron (FCI) limit. Perhaps even more importantly, 
the calculations of [ 2 ^ neglected relativistic, spin-orbit, 
and non-Born-Oppenheimer effects, so accurately esti¬ 
mating Uab initio(i) might Seem impossible. 

However, in my recent benchmark paper [T^, it was 
shown that none of the vibrational energies associated 
with the ab initio potential from [ 2 ^ for the 6(l^nii) state 
of Li 2 , deviated from the empirical potential’s vibrational 
energies by more than 12 cm“^. Since that ab initio po¬ 
tential used the same basis set and method as their po¬ 
tentials for the electronic states approaching 2S + 3P, I 
aimed to make I4rLR(*) — Vab initio(i) less than 15 cm“^ 
for all i except at very small internuclear distances near 
the r = 0 singularity where the inner wall of the potential 
rapidly increases, crosses the dissociation limit, and then 
attains extremely large energy values. The exact values 
used for Uab iniUoii) that were used are presented in Ta¬ 
bles Hvl ED Em] and B Furthermore, there are places 
in which it was desirable to make Uab initio{i) smaller 
than 15 cm“^. This was in places where the potentials 
from had features such as tiny second minima, tiny 
shelves, or any other type of abrupt change. The sub¬ 
sections (below) for each electronic state will describe in 
detail the nature of these features and how this affected 
the choice of Uab initio{i) (once again the exact values are 
given in Tables |T3 and |X|. 

The MLR model was fitted to the points from [ 2 ^ with 
various manually adjusted values of the MLR parameters 
(iV^,p, g, rref) in search for the lowest dd according to 
Eq. [T7| such that increasing Np no longer reduced dd 
signihcantly further than the best dd obtained with the 
previous increase in Np. Details for each electronic state 
are described in the subsections below which focus on 
each state. 


A. The 36(3'^nu) state 

The first b state of Li 2 dissociates to the 2S + 2P 
asymptote. A very detailed analysis of theory vs exper- 
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Table III. Parameters defining the MLRp'q^(Af/ 3 ) potentials fitted to ab initio points from [^, and with long-range functions u{r) 
defined according to the descriptions in section ITCl and with long-range constants presented in Table HIl Numbers in parentheses 
are 95% confidence limit nncertainties in the last digit (s) shown, calculated from the least-squares fitting procedure. 



3b(3®n„) 

MLR; 

i;?(i7) 



3B(3^n,) 

MLR|;i(4) 


3C{3^Ilg) 

MLR?;g(8) 

3d(3®n9) 

MLR|;i(4) 

Se 

5765.593 cm-i 



De 

5368.8(38) cm-i 

De 

4066.0467 cm-^ 

De 6045.135 cm-i 

Te 

3.984 527(23) A 



Te 

3.165 8(18) A 

re 

3.184 402(23) A 

re 3.137 142(23) A 

00 

0.378369938 

010 

-3213.8074 

00 

-0.8509 

00 

0.78516117 

00 0.78516117 

01 

0.918196351 

011 

1673.5923 

01 

-0.592 

01 

2.3951022 

01 2.39510218 

P 2 

2.08782371 

012 

6482.429 

P 2 

-0.038 

02 

6.950787 

02 6.9507868 

03 

0.96103485 

013 

-1705.759 

03 

0.74 

Ps 

-4.68625 

03 -4.686246 

04 

-46.0405588 

014 

-7329.75 

04 

0.57 

04 

-5.17796 

04 -5.17796 

05 

-117.75187 

015 

822.19 



Pb 

5.22029 


06 

271.27645 

016 

4297.08 



Pg 

7.5253 


Pr 

880.09337 

017 

-122.54 



Pr 

-2.100 



iment for the first b state of Li 2 was recently reported 
which summarized 14 different experiments pro- 
viding new information on the 6 (l^nii) state since 1983 
P, I8fil4^ , and mentioned several other papers that in¬ 
volved this state without providing information on any 
new levels. Due to the spin-orbit coupling between b 
states of alkali dimers and their respective yl(^S+) states, 
recent experimental and theoretic al pa pers ha ve st udied 
the lowest b state of R b-?. 1^ . 1^ . llOOl l , Na Cs [lOl| . KCs 
[Tn^ . RbCs Pol, Cs 2 Poland NaK pf)| . 

It is thus surprising that no experiments have been 
reported for the second b state of Li 2 , which dissociates to 
the 2P + 2P asymptote. The present paper is concerned 
with the third b state, which dissociates to 2S + 3P. The 
ab initio potential for 36(3^nu) from pl has a small 
shelf-like feature before the minimum, and another much 
longer one closer to dissociation (see Fig. [T|). The first 
shelf is located between u = 0 and v = 1, and it lasts 
from about 3.2 — 3.6 A. The second shelf starts after 
V = 27 and lasts from about 6.4 — 8 A. Despite these 
fairly pronounced shelf features, at the resolution of the 
ab initio points (which is about 0.1 A), the 36 state only 
has one minimum! 

In Section El it was mentioned that, with the excep¬ 
tion of points at very small values of r, the goal was to 
match all ab initio points of to within ±15 cm“^ 
(and points at larger values of r even better), since this 
was about the level of accuracy found when comparing 
ab initio points psj ] to an empirical potential for 6 (l^n„) 
in p| . Preliminary fits used such a weighting scheme 
for the least-squares fitting, except with points com¬ 
prising the two shelves mentioned in the previous para¬ 
graph, weighted with much smaller uncertainties. This 
was especially important for the second shelf, which only 
spanned a range of < 100 cm“^: Because if the discrep¬ 
ancy between the MLR and ab initio was —15 cm“^ at 
one point, and ±15 cm“^ at another point, the 30 cm“^ 
range of discrepancy would constitute a significant por¬ 
tion of the range of the entire shelf itself. After these 


preliminary fits, it was found that in order to get the 
MLR matching the original data with the desired preci¬ 
sion, it helped to decrease the uncertainties on the non¬ 
shelf ab initio points to slightly below ±15 cm“^. The 
best fits were found with the weights shown in Table 
IIVI Once these final weights were chosen, fits were per¬ 
formed for 252 different combinations of the MLR pa¬ 
rameters {Np,p,q,rj;ef), with 4 < < 17, 6 < p < 8 , 

2 < 9 < 8 , (5 < Tref < 6.5) A, though not every point 
in the convex hull formed by these ranges was used. For 
example, some fits were done with {p,q) = ( 6 , 2 ) and 
some fits were done with (p, q) = ( 8 , 6 ) but it did not 
seem necessary to do fits with (p, 9 ) = (8,2). The best 
fit was found with (Ag,p, 9 , rj-ef) = (17, 6 , 7, 5.9 A), and 
had dd = 1.548, while the best fit with Ag = 16 was 
with {Np,p, 9 , Tj-ef) = (16,6, 6 , 5.3 A) and had dd = 2.790. 
Apart from in the inner wall of the potential, the worst 
discrepancy between the MLR and the original points 
for this Ag = 16 case it was > 30 cm“^, while for this 
Ag = 17 case was < 11 cm“^ so it was quite easy to 
select the Ag = 17 fit. Since this Ag = 17 fit satisfied all 
of our desiderata, Ag = 18 fits were not explored. 

The final MLR parameters for the chosen case are given 
in Table IIIII The inset of Fig. [T] shows the long-range 
behavior of the MLR potential and the original ab ini¬ 
tio points of in Le Roy space, and compares them 
to the theoretical long-range potential based on Eq. [5] 
and the long-range constants in Table El The agree¬ 
ment is surprisingly excellent, however after about 17.7 A 
(l/r3 R:! 0.0018A ) we see that the original points dip 

below the theoretical curve, which should not happen 
because Cg is attractive (see Section 4.3 of [H, for ex¬ 
ample). The fact that the MLR potential matches the 
theoretical curve in this regard, is yet another advan¬ 
tage of using an MLR to represent the ab initio points. 
Furthermore, while not shown in Fig. [U the theoretical 
long-range curve without damping, matches the damped 
curve shown in the figure to graphical accuracy at least 
in the range 14A < r < 00 (0 < l/r^ = 0.000 36 A ^), so 
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this conclusion about the original points spuriously dip¬ 
ping below the theoretical curve is true whether or not 
long-range damping is considered. 

B. The 3B(3^nu) state 

The first B state has a barrier before dissociating to the 
25 -I- 2P asymptote, and has therefore been the subject 
of many empirical studies [iMImiL I t has also been 
used to study other states, such as in |94lll2l4ll,il |. The 
second B state potential (sometimes called the 
state” rather than the 2B state) dissociates to 2P + 2P 
and hugs the inside of the first minimum of the 2 A( 1 ^I]+) 
potential, and due to the perturbations between these 2B 
and 2A states, ther e have b een many empirical studies 
of the 2B state l9i,[nMi3. 

The present paper is concerned with third B state. The 
3B state dissociates to 2S + 3P and there has only been 
one experiment which studied t he t hird B state (some¬ 
times called the state”) [la , which was over 55 

years ago! The authors of that work mentioned in their 
paper that they were not able to confidently assign vibra¬ 
tional quantum numbers to their data, and therefore they 
were only able to conclude that Tg < 34140 cm“^ and 
LOg ~ 250 cm“^. In that study, the anharmonic values 
XgUg, equilibrium rotational constants Bg, and dissocia¬ 
tion energies Dg were determined for the 2B state of Li 2 
and for the 2B and 2>B states of Na 2 , but not for the 35 
state of Li 2 . 

It is no surprise that none of the experiments on the 
2B state showed any indication of a barrier in the poten¬ 
tial, because the leading long-range term (Cs) is attrac¬ 
tive [i^. However, it is perhaps surprising that the ab 
initio potential of for the 35 state does not have a 
barrier (at least before the calculations stopped at about 
21 A), because the leading long-range interaction term 
(C 3 ) for this state is repulsive [^. There is however a 
good theoretical explanation for the lack of barrier in the 
35 state. The long-range potential for 35 is not just a 
simple sum of inverse powers (—Ca/r^ — Ce/r® • • •), but 
it is the middle eigenvalue of the 3x3 spin-orbit interac¬ 
tion matrix of Eq. [TOl which involves the 6 a( 6 ^Eii) and 
36(3^nii) states. For large internuclear distances, the 6 a 
state pushes down on the 35 state enough to remove the 
barrier. Indeed, numerical calculations of the eigenvalues 
of the 3x3 matrix show that the potential is attractive at 
all distances (beyond the repulsive inner wall for r <C Vg). 
This 3x3 matrix for 35 is the exact same as the one 
for the first 5 state, which does have a barrier, but for 
the 35 state the C 3 is three orders of magnitude smaller 
than for the first 5 state, and the Cq is one order of 
magnitude bigger than for the first 5 state [i^. This 
highlights the importance of using the 3x3 coupling ma¬ 
trix because a simple inverse power sum with a negative 
leading long-range term (such as the negative C 3 in the 
present case) is guaranteed to have a barrier, and there 
was even a barrier when using the middle eigenvalue of 


the 3x3 coupling matrix for the 2S + 25 values of Cm, 
but the specific values of Cm at 2S + 35 seem to be past 
a bifurcation point at which the barrier is lost. 

Since alkali parent states of 5(^n„) symmetry only 
have one spin-orbit daughter state {^u/g = lu): we do 
not need to worry about defining the MLR long-range 
function u{r) in a piece-wise manner, and it is simply 
defined as the middle eigenvalue of the 3x3 interac¬ 
tion matrix of Eq. [TOl Also, since there are no features 
such as barriers, multiple minima or shelves, the weight¬ 
ing strategy was straightforward. Preliminary fits were 
done with all points from being weighted with un¬ 
certainty ±15 cm“^ if V (r) < —100 cm“^, with ±5 cm“^ 
if V{r) < —20 cm“^, with ±1 cm“^ if V{r) < 1 cm“^ 
and ±0.5 cm“^ for the one point for which —lcm“^ < 
V{r) <0 cm“^. It was then found that most uncertain¬ 
ties could be even further reduced without making the 
fitting too difficult, and that for very small values of r 
it was very difficult to achieve ±15 cm“^ agreement in 
the fit. In the end, the points from [ 2 ^ for the small¬ 
est values of r (with V{r) > —4000cm“^) were weighted 
with uncertainty ±100 cm“^, and all other points were 
weighted with ±I cm“^, except for the last point which 
was weighted with ±0.5 cm“^ because V{r) itself was 
only Ri —0.7 cm“^. These final weights are shown in 
Table ED 

Once these final weights were chosen, fits were per¬ 
formed for 150 different combinations of the MLR pa¬ 
rameters {Ni 3 ,p,q,r^ei), with 3 < Np < 5, 6 < p < 7, 
2 < <7 < 7, (5 < rref < 10 ) A, though not every point in 
the convex hull formed by these ranges was used. The 
best fit with Np = 4 was with (p, g, Cref) = ( 6 , 6 , 6.2 A) 
which had dd = 4.307 while the best fit found using 
Np = 5 was only marginally better {dd = 4.128) and 
no fits were found with Np = 3 that had dd < 10. There¬ 
fore, the choice of MLR model for this electronic state 
was easy to make, and the final parameters are listed in 
Table in 


C. The 3C(3^n9) state 

The first (^(^ng) state dissoci ates to 2S ± 25 and was 
not studied in detail until 1990 |l20l |. This was 11 years 
after the second state of C'(^ng) symmetry (sometimes 
calle d the “G^Hg state” since it wa s giv en this name in 
[l2l| l was studied in detail in 1979 [l22| | . This 2G state 
dissociates to 25 ± 25 and was studied again in a serie s 
of follow-up papers by Bernheim et al. 0, [m [ 113 . 
Impressively, empirical spectroscopic constants have also 
been reported for all R ydbe rg states in the series ndir^Ilg 
for n = 3 - 15 (!) 0.fl^. In the same paper, it was 
determined that the n = 3 state in this series is in fact 
the 2G state. 

The third C state is the subject of the present work, 
since it dissociates to 2S ± 35. It only has one spin- 
orbit daughter state, which has Ig symmetry and cou¬ 
ples to the Ig daughters of the 3c(3^E+) and 3d(3^ng) 
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states. The potential energy of the 3C state’s Ig daugh¬ 
ter is given by the middle eigenvalue of the appropriate 
3x3 matrix. Preliminary fits were done with the same 
weighting scheme as for the 2>B state, except with the 
points surrounding the second minimum weighted more 
strongly (with ±0.5 cm“^), since the depth of this well is 
< 10 cm“^ and therefore it would not be satisfactory to 
fit to these points with an agreement of only ±15 cm“^! 
It was then found that many points could be weighted 
more strongly without making the fitting too difficult. 
The weights were adjusted to ±10 cm“^ for points at the 
smallest values of r {V{r) > —3000 cm“^), to ±5 cm“^ 
for the rest of the points with V{r) < —1000 cm“^, to 
±1 cm“^ for points with V{r) < —1 cm“^, except for 
the points near the second minimum, and the very last 
point for which V{r) —0.75 cm“^. The final weights 
are presented in Table IVIIII 


With these final weights, fits were performed for 
302 different combinations of the MLR parameters 
(7V;3,P, g, r-ref), with 2 < iV ;3 < 11, 6 < p < 11, 
(4 < Tref < 12) A, though not every point in the con¬ 
vex hull formed by these ranges was used. The best fit 
with Np = 8 was with (p, q, r^ef) = (6, 5, 7.5 A) which had 
dd = 1.612, but this potential did not capture the second 
minimum very well (particularly, it approaches the bar¬ 
rier leading to that minimum with the ab initio point of 
y (7.408 A) = —88.76 cm“^ being represented with a dis¬ 
crepancy of > 5.5 cm“^). The best fit with Np = 8 which 
captured y(7.408A) with a discrepancy of < 2 cm“^ 
was with (p, g, rref) = (6, 8, 7.0 A) which had dd = 3.536. 
None of the Np = 7 models that reproduced V (7.408 A) 
with discrepancy < 2 cm“^ had an overall dd < 4, and 
while Njs = 9 fits were found with discrepancies for this 
point < 2 cm“^ and overall dd as low as 2.593, there 
were no points beyond r = 2.6 A for which the Np = 8 
case with dd = 3.536 misrepresented an original point 
by > 18 cm“^ (the highest discrepancy was 17.37 cm“^ 
at 3.175 A, and among these Np = 9 cases, the low¬ 
est discrepancy for this same point was 12.70 cm“^). 
While deciding not to go beyond Np = 8 was not an 
easy choice, there is not much reason to believe that the 
calculation in for y(3.175A) is so precise that rep¬ 
resenting it more closely by ~ 5 cm“^ is worth adding 
an extra parameter. Here it is mentioned that while the 
comparison against the empirical potential in [T^ for the 
lowest b state showed no discrepancy of > 12 cm“^, that 
paper also noted the surprisingly small effect of Born- 
Oppenheimer breakdown in that system, meaning that 
it is likely that the potentials in for other electronic 
states (especially ones approaching 2S — 3P, which seem 
to interact with each other more than the ones approach¬ 
ing the 2S — 2P state) will be accurate to slightly less 
precision than ±12 cm“^. The final MLR parameters for 
the chosen model are listed in Table nni 


D. The 3d(3^ng) state 


The first d state has a potential energy curve which 
approaches the 2S ± 2P asymptote, but the ab initio 
calculations of [2^ indicate that it has no bound states. 
Therefore, it is no surprise that no bound levels have been 
found in experiments on this s tate, thou gh it was indeed 
involved in some experiments 124l - ll26j| . While the pre¬ 
diction in [ 2 ^ that the first d state has no bound levels 
is likely to be true, it should be noted that ab initio pre¬ 
dictions of this sort are not always reliable. The 1995 ab 
initio study of [ 2 ^ predicted that the 1^S“ state would 
have no bound levels, but the 2014 calculations of 
found there to be a dissociation energy of Se =14 cm“^ 
and an equilibrium harmonic frequency of We = 10 cm“^, 
indicating the existence of at least two bound vibrational 
levels! Likewise, the 2006 ab initio study of found 
the l^Au state to not have any bound levels, but the 
earlier 1995 study of [ 2 ^ and the 2014 study of [ 2 ^ both 
predicted De > 3430 cm“^ and We < 255 cm“^. 

The second d state has been studied extensively, ^ec- 
troscopic measurements for 2d were made in [l|, y, 

H, [^, [ 93 , and 2d was also used i n various other experi¬ 
mental studies such as [9l|, IH, SlUMli. 

The focus of this paper is on the third d state, which 
is the only A = H state dissociating to 2S ± 3P for which 
rovibrationally resolved spectra have been measured, but 
only 13 lines were observed (with v = 6,7, 8,10) 0. Hy- 
perfine structure was also studied experimentally for 3d 
in Q, but it was only for the N = 6,8 levels of u = 8, 
which had already been studied without focus on hyper- 
fine structure in [6|. Finally, the 3d state was involved in 
the experiments of Q, but the focus of that study was 
not the 3d state. 

Since the 3d state has four spin-orbit daughter states 
(Og , 0“, Ig, 2g, analogous with the 36 state), we treat the 
2g symmetry in this paper, since there is only one state 
dissociating to 2S ± 3P with 2g symmetry and therefore 
the long-range potential is a simple sum of inverse-power 
terms rather than a complicated 2 x 2 or 3 x 3 interaction 
matrix. The fitting strategy was very similar to what 
it was for the 3C state, except the 3d state seemed to 
require stronger weighting of the points near the second 
minimum, and weaker weighting of other points. These 
final weights are shown in Table IVII 

Once these final weights were chosen, fits were per¬ 
formed for 238 different combinations of the MLR pa¬ 


rameters {Np,p,q,r„ef), with 3 < < 11, 5 < p < 8, 

2 < g < 10, (5 < Tref < 7.5) A, though not every point in 
the convex hull formed by these ranges was used, and it 
is noted that p must be > 6 in order to ensure the correct 
long-range behavior 0 , but fits with p = 5 were still in¬ 
structive to better understand the model dependence for 
this potential. The best fit with = 4 and p > 6 was 
with (p, g, Tref) = (6, 7, 7.2 A) which had dd = 4.292, only 
0.043 higher than the best fit found with p = 5. Fits with 
Np = 5 had dd values as low as 2.972 but the mentioned 
Np = 4 case did not misrepresent any of the original 
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points beyond 2.434 A) by > 13 cni“^, so all desiderata 
were satisfied without resorting to = 5. No Np = 3 
cases had dd < 10, so the choice of MLR model for this 
electronic state was easy to make, and the final parame¬ 
ters are listed in Table EH 

III. CONCLUSION 

Analytic MLR potentials were htted to the ab initio 
points from and with correct long-range behavior in¬ 
corporated according to effects of spin-orbit coupling de¬ 
scribed in Eqs. 151 [TUI [TT] and the long-range constants in 
Table El Despite the potentials from [2^ having unusual 
features such as multiple minima, barriers, and shelves, 
which have never been described by an MLR-type model 
before, all of these features were successfully captured 
with the MLR model. This answers an age-old ques¬ 
tion of whether or not fully analytic potentials can have 
the flexibility needed in order to capture such features. 
Pashov’s 2008 paper “Pointwise and analytic potentials 
for diatomic molecules. An attempt for critical compari¬ 
son” described lack of flexibility as one of the three 
drawbacks of analytic potentials, and suggested that the 
MLR model may not be able to capture double minima 
or shelf-like features. Five years earlier in 2003, Huang 
and Le Roy suggested in that Pashov’s pointwise ap¬ 
proach would be the method of choice for potentials such 
as those described in this paper: 

“A particular strength of [Pashov’s point- 
wise] model is the fact that it has more lo¬ 
cal flexibility than do fully analytical poten¬ 
tial function forms, in that a shift of one po¬ 
tential point has only a modest effect on the 
function outside its immediate neighborhood. 

This would tend to make [Pashov’s pointwise] 
model the method of choice for cases where 
the potential has substantial local structure 
or undergoes an abrupt change of character 
on a small fraction of the overall interval, such 
as occurs near an avoided curve crossing. In 
contrast, a change in one of the parameters 
defining a [fully analytic] such as our DELR 
function will in general affect the potential 
across the whole domain. This makes the 
parameters defining [fully analytic] potentials 
very highly correlated and can give rise to dif¬ 
ficulty in achieving full unique convergence in 
a fit.” 

At the time when this quote was written, the MLR model 
did not exist yet, and at the time of Pashov’s paper [s^] . 
only a primitive (less flexible) form of the MLR existed, 
which was used in just four simple cases of ground elec¬ 
tronic states (Sa - I^ . It is possible that the notion that 
analytic potentials cannot capture such special features 
may be attributed to the lack of diversity in attempted 
applications at that early stage in time, and the lack of 


some of the newer MLR features which were introduced 
in [T^ and for increasing flexibility, and in for 
correcting the long-range behavior. 

This paper also represents, to the best of my knowl¬ 
edge, the most detailed study of analytic potentials for 
excited electronic molecular states. 
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Table IV. Quality of MLRe'7(17) fit to original ab initio en¬ 
ergies of for the 3fc(3®n2„) state. When an isotopologue 
listed in Table |V] has at least one vibrational energy level be¬ 
yond the range of ab initio data available, approximate dis- 


Table V. Vibrational energies in cm ^ for the 3fe(3®n2„) state 
predicted by the MLRg;7(17) from Table Hill 

V ®’®Li2 ®-"Li2 "'"Li2 “’“Li2 


tances are given for the outer classical turning points of the 


jrresponding vibrational wavefunctions. 


Se 

5 765.593 5 765.593 5 765.593 5 765.593 







0 

-5 654.932 -5658.488 -5662.222 -5661.726 

r 


Original 

Weight 

MLR - Original 

MLR 

1 

-5 493.765 -5500.003 -5506.570 -5505.696 

fAl 






2 

-5 385.710 -5397.683 -5356.570 -5408.260 



[Clll J 

[Clll J 

[Clll J 

[Clll J 

3 

-5 241.269 -5260.427 -5280.073 -5277.492 







4 

-5 083.976 -5110.717 -5138.177 -5134.567 

2.328 

380 

-1 501.30 

40.00 

195.839 

-1305.461 

5 

-4914.852 -4980.180 -4984.997 -4980.318 

2.434 

215 

-2 661.44 

40.00 

9.542 

-2651.898 

6 

-4 737.691 -4780.180 -4823.968 -4818.203 

2.540 

051 

-3 556.46 

40.00 

-79.794 

-3636.254 

7 

-4 554.829 -4605.151 -4657.110 -4650.263 

2.645 

886 

-4 232.00 

40.00 

-104.149 

-4336.149 

8 

-4 368.240 -4426.227 -4486.207 -4478.297 

2.751 

721 

-4 727.14 

40.00 

-89.498 

-4816.638 

9 

-4179.292 -4244.786 -4312.615 -4303.664 

2.857 

557 

-5 075.00 

20.00 

-57.523 

-5132.523 

10 

-3 988.756 -4061.679 -4137.249 -4127.274 

2.963 

392 

-5 304.35 

10.00 

-25.034 

-5329.384 

11 

-3 796.972 -3877.308 -3960.593 -3949.597 

3.069 

228 

-5 440.87 

0.70 

-2.753 

-5443.623 

12 

-3 604.218 -3691.902 -3782.874 -3770.860 

3.175 

063 

-5 508.47 

0.70 

5.647 

-5502.823 

13 

-3 411.033 -3505.823 -3604.333 -3591.316 

3.280 

899 

-5 530.85 

0.70 

3.179 

-5527.671 

14 

-3 218.297 -3319.731 -3404.333 -3411.444 

3.386 

734 

-5 533.92 

0.70 

-1.971 

-5535.891 

15 

-3 027.108 -3134.538 -3246.854 -3231.976 

3.492 

570 

-5 543.14 

0.70 

-2.801 

-5545.941 

16 

-2 838.477 -2951.238 -3069.490 -3053.803 

3.598 

405 

-5 576.94 

0.70 

1.583 

-5575.357 

17 

-2 652.930 -2770.592 -2894.177 -2877.765 

3.704 

240 

-5 634.66 

0.70 

3.546 

-5631.114 

18 

-2 470.302 -2592.818 -2721.446 -2704.363 

3.810 

076 

-5 700.51 

0.70 

0.373 

-5700.137 

19 

-2 290.105 -2417.585 -2551.310 -2533.558 

3.915 

911 

-5 750.55 

0.70 

-2.666 

-5753.216 

20 

-2 112.256 -2244.505 -2383.395 -2364.958 

4.021 

747 

-5 761.08 

0.70 

-0.047 

-5761.127 

21 

-1937.405 -2073.700 -2217.430 -2198.327 

4.127 

582 

-5 708.41 

0.90 

2.634 

-5705.776 

22 

-1766.631 -1905.899 -2053.644 -2033.962 

4.233 

418 

-5 578.48 

0.90 

-0.382 

-5578.862 

23 

-1600.818 -1742.086 -1892.754 -1872.628 

4.339 

253 

-5 376.12 

1.00 

-3.307 

-5379.427 

24 

-1440.124 -1582.955 -1735.630 -1715.198 

4.445 

089 

-5 120.65 

2.00 

1.698 

-5118.952 

25 

-1284.340 -1428.527 -1582.844 -1562.186 

4.550 

924 

-4 832.92 

4.00 

9.143 

-4823.777 

26 

-1 134.583 -1278.624 -1434.376 -1413.503 

4.656 

759 

-4 528.51 

4.00 

7.006 

-4521.504 

27 

-996.778 -1134.368 -1290.038 -1269.027 

4.762 

595 

-4 217.95 

4.00 

-4.138 

-4222.088 

28 

-914.777 -1001.099 -1150.711 -1129.946 

4.868 

430 

-3 908.28 

4.00 

-9.716 

-3917.996 

29 

-882.576 -916.881 -1020.440 -1001.313 

4.974 

266 

-3 603.65 

1.00 

-1.519 

-3605.169 

30 

-845.026 -887.471 -950.440 -918.124 

5.080 

101 

-3 307.57 

3.00 

10.262 

-3297.308 

31 

-800.327 -851.265 -898.383 -890.415 

5.185 

937 

-3 022.26 

3.00 

9.62 

-3012.64 

32 

-753.924 -830.000 -861.863 -813.832 

5.291 

772 

-2 749.45 

1.00 

-2.696 

-2752.146 

33 

-705.900 -808.055 -820.356 -813.832 

5.609 

278 

-2 018.82 

1.00 

1.364 

-2017.456 

34 

-654.935 -718.573 -780.226 -772.047 

5.820 

949 

-1 611.25 

1.00 

-2.927 

-1614.177 

35 

-602.077 -669.908 -737.158 -728.415 

6.138 

456 

-1 127.31 

1.00 

4.595 

-1122.715 

36 

-548.234 -619.419 -690.970 -681.806 

6.350 

127 

-958.32 

0.30 

-0.912 

-959.232 

37 

-493.706 -567.885 -643.224 -633.418 

6.667 

633 

-935.05 

0.30 

0.886 

-934.164 

38 

-438.926 -515.524 -594.281 -583.945 

6.879 

304 

-931.76 

0.30 

-1.386 

-933.146 

39 

-384.474 -462.746 -544.287 -533.563 

7.196 

810 

-925.40 

0.30 

1.252 

-924.148 

40 

-330.903 -410.071 -493.717 -482.653 

7.408 

481 

-917.28 

0.30 

-1.012 

-918.292 

41 

-278.814 -357.964 -443.010 -431.682 

7.725 

987 

-895.77 

1.50 

3.616 

-892.154 

42 

-228.913 -306.931 -392.539 -381.054 

7.937 

658 

-872.50 

1.50 

5.113 

-867.387 

43 

-182.018 -257.565 -342.733 -331.205 

8.255 

164 

-817.85 

1.50 

-2.201 

-820.051 

44 

-139.060 -210.538 -362.733 -282.642 

8.466 

835 

-764.30 

4.00 

-6.928 

-771.228 

45 

-101.066 -166.620 -247.119 -235.929 

8.784 

342 

-656.98 

4.00 

-5.076 

-662.056 

46 

-69.047 -126.682 -202.462 -191.711 

8.996 

013 

-572.26 

4.00 

9.21 

-563.05 

47 

-43.749 -91.656 -160.811 -150.725 

9.260 

601 

-464.06 

4.00 

8.64 

-455.42 

48 

-25.299 -62.409 -122.952 -113.784 

9.525 

190 

-364.42 

4.00 

8.00 

-356.42 

49 

-13.022 -39.497 -89.714 -81.728 

9.789 

778 

-280.36 

4.00 

7.29 

-273.07 

50 

-5.669 -22.881 -61.855 -55.277 

10.054 

367 

-213.64 

4.00 

6.54 

-207.1 

51 

-1.855 -11.834 -39.854 -34.795 

10.318 

956 

-163.16 

4.00 

5.78 

-157.38 

52 

-0.324 -5.193 -23.675 -20.079 

10.583 

544 

-125.85 

2.00 

5.04 

-120.81 

53 

-0.002 -1.723 -12.691 -10.348 

11.112 

721 

-78.44 

2.00 

4.35 

-74.09 

54 

-0.312 -5.885 -4.517 

11.641 

899 

-52.11 

2.00 

3.69 

-48.42 

55 

-0.003 -3.885 -1.482 

12.171 

076 

-36.52 

2.00 

3.16 

-33.36 

56 

-0.493 -0.261 

12.700 

253 

-26.43 

2.00 

2.72 

-23.71 

57 

-0.002 

13.229 

430 

-19.84 

2.00 

2.37 

-17.47 



13.758 

607 

-15.24 

5.50 

2.09 

-13.15 



14.287 

785 

-11.94 

5.50 

1.88 

-10.06 



14.816 

962 

-9.53 

5.50 

1.71 

-7.82 



15.346 

139 

-7.55 

5.00 

1.53 

-6.02 



15.875 

316 

-6.24 

5.00 

1.37 

-4.87 



16.933 

671 

-4.04 

5.00 

1.21 

-2.83 



17.992 

025 

-2.73 

5.00 

1.08 

-1.65 



19.050 

380 

-1.85 

5.00 

0.96 

-0.89 



20.108 

734 

-1.19 

5.00 

0.86 

-0.33 



21.167 

088 

-0.75 

0.50 

0.73 

-0.02 



65. 


— 

— 

— 

-0.003 



70. 


— 

— 

— 

-0.002 
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Table VI. Quality of MLR6;6(4) fit to original ab initio ener¬ 
gies of [ 2 ^ for the 3B(3^ni„). When an isotopologue listed 
in Table En has at least one vibrational energy level beyond 
the range of ab initio data available, approximate distances 
are given for the outer classical turning points of the corre¬ 
sponding vibrational wavefunctions. 


r Original 

[A] [cm-i] 

Weight 

[cm-i] 

MLR - Original 
[cm-i] 

MLR 

[cm-i] 

2.223 -305.82 

100 

798.584 

492.764 

2.328 -1574.16 

100 

354.615 

-1219.545 

2.434 -2650.69 

100 

133.127 

-2517.563 

2.540 -3531.22 

100 

45.218 

-3486.002 

2.646 -4198.64 

1.00 

5.965 

-4192.675 

2.752 -4684.78 

1.00 

-7.168 

-4691.948 

2.858 -5020.35 

1.00 

-7.054 

-5027.404 

2.963 -5232.81 

1.00 

-1.234 

-5234.044 

3.069 -5344.96 

1.00 

5.069 

-5339.891 

3.175 -5375.90 

1.00 

8.448 

-5367.452 

3.281 -5342.54 

1.00 

7.751 

-5334.789 

3.387 -5259.14 

1.00 

2.771 

-5256.369 

3.493 -5140.19 

1.00 

-3.585 

-5143.775 

3.598 -4997.31 

1.00 

-8.898 

-5006.208 

3.704 -4841.04 

1.00 

-9.908 

-4850.948 

3.810 -4677.75 

1.00 

-5.925 

-4683.675 

3.916 -4509.20 

1.00 

0.439 

-4508.761 

4.022 -4335.59 

1.00 

6.090 

-4329.500 

4.128 -4157.16 

1.00 

8.759 

-4148.401 

4.233 -3975.87 

1.00 

8.547 

-3967.323 

4.339 -3793.49 

1.00 

5.762 

-3787.728 

4.445 -3612.64 

1.00 

1.858 

-3610.782 

4.551 -3435.53 

1.00 

-2.007 

-3437.537 

4.657 -3264.12 

1.00 

-4.836 

-3268.956 

4.763 -3099.95 

1.00 

-6.043 

-3105.993 

4.868 -2943.90 

1.00 

-5.667 

-2949.567 

4.974 -2796.86 

1.00 

-3.664 

-2800.524 

5.080 -2658.59 

1.00 

-0.931 

-2659.521 

5.186 -2529.32 

1.00 

2.311 

-2527.009 

5.292 -2407.73 

1.00 

4.610 

-2403.120 

5.609 -2085.10 

1.00 

5.795 

-2079.305 

5.821 -1895.91 

1.00 

1.938 

-1893.972 

6.138 -1632.98 

1.00 

-3.702 

-1636.682 

6.350 -1462.67 

1.00 

-4.605 

-1467.275 

6.668 -1209.84 

1.00 

-1.783 

-1211.623 

6.879 -1046.55 

1.00 

0.657 

-1045.893 

7.197 -822.68 

1.00 

3.489 

-819.191 

7.408 -692.97 

1.00 

3.915 

-689.055 

7.726 -532.10 

1.00 

3.269 

-528.831 

7.938 -445.62 

1.00 

1.783 

-443.837 

8.255 -343.13 

1.00 

-0.567 

-343.697 

8.467 -289.80 

1.00 

-1.852 

-291.652 

8.784 -227.03 

1.00 

-3.294 

-230.324 

8.996 -194.11 

1.00 

-3.988 

-198.098 

9.261 -160.97 

1.00 

-4.249 

-165.219 

9.525 -134.41 

1.00 

-4.343 

-138.753 

9.790 -113.12 

1.00 

-4.130 

-117.250 

10.054 -96.00 

1.00 

-3.622 

-99.622 

10.319 -81.96 

1.00 

-3.106 

-85.066 

10.584 -70.10 

1.00 

-2.845 

-72.945 

11.113 -52.33 

1.00 

-1.958 

-54.288 

11.642 -39.60 

1.00 

-1.360 

-40.960 

12.171 -30.38 

1.00 

-0.904 

-31.284 

12.700 -23.58 

1.00 

-0.585 

-24.165 

13.229 -18.53 

1.00 

-0.326 

-18.856 

13.759 -14.58 

1.00 

-0.287 

-14.867 

14.288 -11.72 

1.00 

-0.128 

-11.848 

14.817 -9.31 

1.00 

-0.232 

-9.542 

15.346 -7.55 

1.00 

-0.189 

-7.739 

15.875 -6.24 

1.00 

-0.093 

-6.333 

16.934 -4.04 

1.00 

-0.275 

-4.315 

17.992 -2.73 

1.00 

-0.292 

-3.022 

19.050 -1.85 

1.00 

-0.310 

-2.160 

20.109 -1.19 

1.00 

-0.385 

-1.575 

21.167 -0.75 

0.50 

-0.421 

-1.171 


Table VII. Vibrational energies in cm ^ for the 3i?(3^ni„) 
state predicted by the MLRg;6(4) from Table Hill 

V ®’«Li2 «'^Li2 ^’^Li2 “’llLi2 


De 5368.8 5368.8 5368.8 5368.8 

0 -5259.702 -5263.694 -5267.840 -5267.292 

1 -5047.141 -5058.712 -5070.746 -5069.155 

2 -4838.570 -4857.433 -4877.067 -4874.470 

3 -4634.049 -4659.912 -4686.855 -4683.290 

4 -4433.598 -4466.170 -4500.132 -4495.636 

5 -4237.212 -4276.205 -4316.896 -4311.508 

6 -4044.874 -4090.002 -4137.136 -4130.892 

7 -3856.578 -3907.553 -3960.843 -3953.781 

8 -3672.341 -3728.865 -3788.017 -3780.174 

9 -3492.217 -3553.977 -3618.683 -3610.100 

10 -3316.304 -3382.962 -3452.895 -3443.613 

11 -3144.757 -3215.941 -3290.743 -3280.808 

12 -2977.780 -3053.079 -3132.356 -3121.818 

13 -2815.625 -2894.584 -2977.902 -2966.817 

14 -2658.578 -2740.700 -2827.582 -2816.010 

15 -2506.934 -2591.690 -2681.624 -2669.630 

16 -2360.954 -2447.812 -2540.262 -2527.916 

17 -2220.808 -2309.273 -2403.714 -2391.085 

18 -2086.496 -2176.175 -2272.140 -2259.292 

19 -1957.783 -2048.452 -2145.593 -2132.577 

20 -1834.183 -1925.822 -2023.966 -2010.813 

21 -1715.028 -1807.793 -1906.966 -1893.682 

22 -1599.606 -1693.743 -1794.128 -1780.695 

23 -1487.302 -1583.044 -1684.891 -1671.279 

24 -1377.678 -1475.167 -1578.702 -1564.878 

25 -1270.487 -1369.744 -1475.101 -1461.040 

26 -1165.652 -1266.570 -1373.763 -1359.458 

27 -1063.233 -1165.583 -1274.505 -1259.962 


28 

-963.397 

-1066.840 

-1177.266 

-1162.507 

29 

-866.401 

-970.488 

-1082.087 

-1067.149 

30 

-772.574 

-876.755 

-989.090 

-974.022 

31 

-682.312 

-785.929 

-898.465 

-883.330 

32 

-596.065 

-698.356 

-810.456 

-795.329 

33 

-514.332 

-614.434 

-725.359 

-710.329 

34 

-437.645 

-534.599 

-643.512 

-628.681 

35 

-366.550 

-459.324 

-565.294 

-550.774 

36 

-301.572 

-389.096 

-491.112 

-477.030 

37 

-243.166 

-324.395 

-421.395 

-407.885 

38 

-191.657 

-265.654 

-356.577 

-343.778 

39 

-147.184 

-213.218 

-297.065 

-285.115 

40 

-109.665 

-167.284 

-243.214 

-232.239 

41 

-78.804 

-127.871 

-195.278 

-185.379 

42 

-54.137 

-94.805 

-153.373 

-144.621 

43 

-35.092 

-67.741 

-117.455 

-109.880 

44 

-21.035 

-46.216 

-87.316 

-80.910 

45 

-11.285 

-29.692 

-62.618 

-57.343 

46 

-5.113 

-17.586 

-42.931 

-38.724 

47 

-1.738 

-9.275 

-27.767 

-24.548 

48 

-0.325 

-4.092 

-16.603 

-14.276 

49 

-0.007 

-1.326 

-8.884 

-8.884 

50 


-0.222 

-4.015 

-4.015 

51 


-0.003 

-1.362 

-7.331 

52 



-0.256 

-3.100 

53 



-0.007 

-0.929 
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Table VIII. Quality of MLRg'g(8) fit to original ab initio en¬ 
ergies of [2^ for the 8(7(3^1113) state. 


r Original 

[A] [cm~i] 

Weight 

[cm-i] 

MLR - Original 
[cm-i] 

MLR 

[cm-i] 

2.328 -331.50 

10.00 

35.693 

-295.807 

2.434 -1429.53 

10.00 

-7.764 

-1437.294 

2.540 -2276.70 

10.00 

-22.184 

-2298.884 

2.645 -2916.47 

10.00 

-20.557 

-2937.027 

2.751 -3385.05 

5.00 

-11.667 

-3396.717 

2.857 -3712.29 

5.00 

-1.457 

-3713.747 

2.963 -3915.52 

5.00 

-0.992 

-3916.512 

3.069 -4031.62 

5.00 

4.084 

-4027.536 

3.175 -4078.81 

5.00 

14.101 

-4064.709 

3.280 -4049.18 

5.00 

6.865 

-4042.315 

3.386 -3981.36 

5.00 

9.450 

-3971.910 

3.492 -3867.67 

5.00 

4.697 

-3862.973 

3.598 -3723.48 

5.00 

0.049 

-3723.431 

3.704 -3556.02 

5.00 

-3.998 

-3560.018 

3.810 -3371.66 

5.00 

-6.935 

-3378.595 

3.915 -3176.11 

5.00 

-8.227 

-3184.337 

4.021 -2973.97 

5.00 

-7.885 

-2981.855 

4.127 -2768.98 

5.00 

-6.346 

-2775.326 

4.233 -2564.65 

5.00 

-3.854 

-2568.504 

4.339 -2364.05 

5.00 

-0.696 

-2364.746 

4.445 -2169.38 

5.00 

2.391 

-2166.989 

4.550 -1982.83 

5.00 

5.140 

-1977.690 

4.656 -1805.27 

5.00 

6.500 

-1798.770 

4.762 -1638.69 

5.00 

7.099 

-1631.591 

4.868 -1483.30 

5.00 

6.415 

-1476.885 

4.974 -1339.76 

5.00 

4.958 

-1334.802 

5.080 -1208.08 

5.00 

3.099 

-1204.981 

5.185 -1087.37 

5.00 

0.726 

-1086.644 

5.291 -977.63 

1.00 

-1.168 

-978.798 

5.609 -707.90 

1.00 

-0.784 

-708.684 

5.820 -568.75 

1.00 

1.689 

-567.061 

6.138 -405.90 

1.00 

0.725 

-405.175 

6.350 -320.52 

1.00 

-2.067 

-322.587 

6.667 -218.69 

1.00 

-0.594 

-219.284 

6.879 -166.01 

1.00 

2.617 

-163.393 

7.196 -109.83 

1.00 

1.706 

-108.124 

7.408 -88.76 

0.50 

-1.103 

-89.863 

7.725 -80.64 

0.50 

-0.478 

-81.118 

7.937 -82.83 

0.50 

0.689 

-82.141 

8.255 -87.00 

0.50 

0.527 

-86.473 

8.466 -88.32 

0.50 

-0.061 

-88.381 

8.784 -87.44 

0.50 

-0.480 

-87.920 

8.996 -85.25 

1.00 

-0.356 

-85.606 

9.260 -81.08 

1.00 

-0.110 

-81.190 

9.525 -76.03 

1.00 

0.225 

-75.805 

9.789 -70.32 

1.00 

0.334 

-69.986 

10.054 -64.40 

1.00 

0.321 

-64.079 

10.318 -58.47 

1.00 

0.196 

-58.274 

10.583 -52.98 

1.00 

0.266 

-52.714 

11.112 -42.67 

1.00 

0.025 

-42.645 

11.641 -34.11 

1.00 

-0.050 

-34.160 

12.171 -27.09 

1.00 

-0.159 

-27.249 

12.700 -21.60 

1.00 

-0.120 

-21.720 

13.229 -17.21 

1.00 

-0.140 

-17.350 

13.758 -13.92 

1.00 

-0.003 

-13.923 

14.287 -11.28 

1.00 

0.034 

-11.246 

14.816 -9.09 

1.00 

-0.061 

-9.151 

15.346 -7.33 

1.00 

-0.151 

-7.481 

15.875 -6.02 

1.00 

-0.140 

-6.160 

16.933 -4.04 

1.00 

-0.194 

-4.234 

17.992 -2.73 

1.00 

-0.253 

-2.983 

19.050 -1.85 

1.00 

-0.290 

-2.140 

20.108 -0.75 

0.50 

-0.415 

-1.165 


Table IX. Vibrational energies in cm ^ for the 8(7(3^1113) 
state predicted by the MLRg;g(8) from Table Hill 

V ®'«Li2 ^'^Li2 “'iiLb 


De 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


-3957.890 

-3745.119 

-3445.119 

-3324.020 

-3116.570 

-2911.754 

-2710.006 

-2511.771 

-2317.519 

-2127.744 

-1942.968 

-1763.742 

-1590.632 

-1424.202 

-1264.977 

-1113.398 

-969.785 

-834.340 

-707.248 

-588.846 

-479.738 

-380.571 

-291.279 

-210.843 

-179.563 

-139.563 

-75.962 

-59.183 

-44.917 

-32.408 

-23.147 

-13.147 

-6.897 


-3961.862 

-3761.862 

-3552.904 

-3350.648 

-3150.393 

-2952.521 

-2757.416 

-2565.475 

-2377.112 

-2192.761 

-2092.761 

-1837.961 

-1668.499 

-1505.005 

-1347.972 

-1197.836 

-1054.947 

-919.540 

-791.774 

-671.841 

-560.115 

-457.196 

-363.593 

-279.026 

-202.533 

-134.770 

-95.896 

-75.896 

-59.831 

-45.935 

-33.705 

-23.705 

-14.533 

-8.064 

-3.706 


-3965.987 

-3765.987 

-3572.894 

-3378.354 

-3185.613 

-2995.010 

-2806.885 

-2606.885 

-2439.475 

-2260.933 

-2086.359 

-1916.175 

-1750.820 

-1590.747 

-1436.402 

-1288.205 

-1188.205 

-1011.602 

-883.647 

-762.784 

-649.211 

-543.311 

-445.642 

-356.603 

-275.840 

-202.474 

-137.129 

-107.129 

-77.432 

-57.432 

-48.506 

-36.508 

-25.881 

-17.012 

- 10.101 

-5.173 


-3965.442 

-3765.442 

-3570.251 

-3374.690 

-3180.954 

-2989.387 

-2800.335 

-2600.335 

-2431.207 

-2251.888 

-2076.600 

-1905.775 

-1739.859 

-1579.313 

-1424.589 

-1276.109 

-1134.234 

-934.234 

-871.271 

-750.496 

-637.121 

-531.558 

-434.379 

-345.930 

-265.738 

-192.963 

-128.680 

-108.680 

-75.195 

-60.065 

-46.489 

-36.489 

-24.164 

-15.608 

-9.034 

-4.447 
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Table X. Quality of MLRg'9(4) fit to original ab initio energies 
of [2^ for the 3d(3^n2g) state. When an isotopologue listed 
in Table lU has at least one vibrational energy level beyond 
the range of ab initio data available, approximate distances 
are given for the outer classical turning points of the corre¬ 
sponding vibrational wavefunctions. 

r Original Weight MLR - Original MLR 
[A] [cm“^] [cm“^] [cin“^] 


Table XL Vibrational energies in cm ^ for the 3d(3^n2g) state 
predicted by the MLRg;9(4) from Table Hill 


V 6'®Li2 

®’^Li2 

^’^Li2 

^l’“Li2 

De 6 045.135 

6 045.135 

6 045.135 

6 045.135 

0 -5909.082 

-5914.058 

-5919.226 

-5918.543 

1 -5641.796 

-5656.375 

-5671.532 

-5669.528 

2 -5377.524 

-5401.482 

-5371.532 

-5369.528 

3 -5116.452 

-5149.547 

-5183.994 

-5179.438 

4 -4858.751 

-4949.547 

-4944.436 

-4938.652 

5 -4604.574 

-4655.149 

-4707.858 

-4700.882 

6 -4354.070 

-4412.959 

-4474.377 

-4466.246 

7 -4107.378 

-4174.278 

-4244.107 

-4234.860 

8 -3864.636 

-3939.232 

-4017.158 

-4006.835 

9 -3625.987 

-3707.944 

-3793.641 

-3782.284 

10 -3391.579 

-3480.547 

-3573.669 

-3561.322 

11 -3161.574 

-3257.179 

-3357.363 

-3344.073 

12 -2936.152 

-3037.995 

-3144.853 

-3130.670 

13 -2715.518 

-2823.167 

-2936.284 

-2921.261 

14 -2499.912 

-2612.894 

-2731.819 

-2716.013 

15 -2289.614 

-2407.408 

-2531.648 

-2515.121 

16 -2084.961 

-2206.980 

-2331.648 

-2318.811 

17 -1886.358 

-2011.936 

-2145.098 

-2118.811 

18 -1694.302 

-1822.668 

-1959.285 

-1941.046 

19 -1509.404 

-1639.650 

-1778.914 

-1760.288 

20 -1332.434 

-1463.471 

-1604.431 

-1585.534 

21 -1164.385 

-1294.865 

-1436.381 

-1417.351 

22 -1006.574 

-1134.778 

-1275.444 

-1256.450 

23 -860.858 

-984.465 

-1122.490 

-1103.740 

24 -730.125 

-845.692 

-978.663 

-960.432 

25 -620.025 

-721.209 

-845.563 

-828.230 

26 -600.025 

-616.435 

-725.642 

-709.810 

27 -536.279 

-586.435 

-623.550 

-610.540 

28 -502.410 

-537.243 

-582.928 

-576.909 

29 -464.979 

-505.003 

-542.928 

-536.909 

30 -425.955 

-469.139 

-511.362 

-506.341 

31 -380.025 

-431.738 

-477.652 

-476.341 

32 -344.186 

-411.738 

-442.209 

-435.772 

33 -302.862 

-353.180 

-405.126 

-398.323 

34 -262.005 

-313.325 

-367.188 

-360.082 

35 -222.176 

-273.774 

-328.906 

-321.577 

36 -183.954 

-235.025 

-290.723 

-283.260 

37 -163.228 

-197.591 

-253.083 

-245.580 

38 -113.198 

-162.009 

-216.432 

-208.993 

39 -104.815 

-128.846 

-181.237 

-173.977 

40 -79.439 

-98.692 

-117.188 

-150.689 

41 -39.669 

-69.707 

-107.188 

-110.689 

42 -23.171 

-49.707 

-84.675 

-90.404 

43 -12.113 

-31.768 

-64.675 

-70.404 

44 -5.236 

-18.400 

-44.675 

-40.404 

45 -1.618 

-9.302 

-28.485 

-25.161 

46 

-3.816 

-16.488 

-20.161 

47 

-1.061 

-8.348 

-6.813 

48 


-3.439 

-2.591 

49 


-0.964 

-0.608 


2.222544 -675.63 

150.00 

242.061 

-433.569 

2.328380 -2176.84 

15.00 

96.978 

-2079.862 

2.434215 -3353.22 

10.00 

27.796 

-3325.424 

2.540051 -4256.36 

1.00 

1.853 

-4254.507 

2.645886 -4929.71 

1.00 

-3.720 

-4933.430 

2.751721 -5412.99 

1.00 

-1.158 

-5414.148 

2.857557 -5739.79 

1.00 

2.660 

-5737.130 

2.963392 -5939.29 

1.00 

5.574 

-5933.716 

3.069228 -6034.33 

1.00 

6.128 

-6028.202 

3.175063 -6043.98 

1.00 

4.493 

-6039.487 

3.280899 -5986.48 

1.00 

3.904 

-5982.576 

3.386734 -5867.31 

1.00 

-2.283 

-5869.593 

3.492570 -5704.90 

1.00 

-5.728 

-5710.628 

3.598405 -5506.27 

1.00 

-8.024 

-5514.294 

3.704240 -5278.46 

1.00 

-9.647 

-5288.107 

3.810076 -5028.69 

1.00 

-9.954 

-5038.644 

3.915911 -4762.69 

10.00 

-9.071 

-4771.761 

4.021747 -4485.27 

10.00 

-7.326 

-4492.596 

4.127582 -4200.62 

10.00 

-5.101 

-4205.721 

4.233418 -3912.45 

10.00 

-2.679 

-3915.129 

4.339253 -3624.06 

10.00 

-0.318 

-3624.378 

4.445089 -3338.52 

1.00 

1.889 

-3336.631 

4.550924 -3058.47 

1.00 

3.726 

-3054.744 

4.656759 -2786.32 

1.00 

5.015 

-2781.305 

4.762595 -2524.27 

1.00 

5.575 

-2518.695 

4.868430 -2274.73 

1.00 

5.652 

-2269.078 

4.974266 -2039.01 

1.00 

4.615 

-2034.395 

5.080101 -1819.10 

1.00 

2.755 

-1816.345 

5.185937 -1616.52 

1.00 

0.247 

-1616.273 

5.291772 -1432.60 

1.00 

-2.570 

-1435.17 

5.609278 -997.82 

1.00 

-10.897 

-1008.717 

5.820949 -804.90 

1.00 

-12.271 

-817.171 

6.138456 -638.54 

1.00 

-6.685 

-645.225 

6.350127 -587.84 

0.10 

-1.544 

-589.384 

6.667633 -563.92 

0.10 

1.767 

-562.153 

6.879304 -564.36 

0.10 

0.378 

-563.982 

7.196810 -570.07 

0.10 

-1.508 

-571.578 

7.408481 -570.51 

0.10 

-0.877 

-571.387 

7.725987 -559.75 

0.10 

1.099 

-558.651 

7.937658 -543.07 

1.00 

1.662 

-541.408 

8.466835 -466.48 

1.00 

0.547 

-465.933 

8.784342 -401.07 

1.00 

-1.351 

-402.421 

8.996013 -353.88 

1.00 

-3.111 

-356.991 

9.260601 -296.16 

1.00 

-5.093 

-301.253 

9.525190 -243.05 

1.00 

-7.119 

-250.169 

9.789778 -197.62 

1.00 

-8.143 

-205.763 

10.054367 -160.09 

1.00 

-8.449 

-168.539 

10.318956 -129.80 

1.00 

-8.245 

-138.045 

10.583544 -105.88 

1.00 

-7.519 

-113.399 

11.112721 -72.08 

1.00 

-5.612 

-77.692 

11.641899 -50.57 

1.00 

-3.992 

-54.562 

12.171076 -36.52 

1.00 

-2.767 

-39.287 

12.700253 -26.87 

1.00 

-2.082 

-28.952 

13.229430 -20.50 

1.00 

-1.257 

-21.757 

13.758607 -15.67 

1.00 

-0.966 

-16.636 

14.287785 -12.38 

1.00 

-0.534 

-12.914 

14.816962 -9.75 

1.00 

-0.406 

-10.156 

15.346139 -7.77 

1.00 

-0.300 

-8.070 

15.875316 -6.24 

1.00 

-0.245 

-6.485 

16.933671 -4.04 

1.00 

-0.243 

-4.283 

17.992025 -2.73 

1.00 

-0.189 

-2.919 

19.050380 -1.85 

1.00 

-0.181 

-2.031 

20.108734 -1.19 

1.00 

-0.251 

-1.441 

21.167088 -0.75 

0.50 

-0.290 

-1.040 

21.5 

— 

— 

-0.964 

23. 

— 

— 

-0.608 
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